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I. INTRODUCTION 

The use of radioactive nuclear beams produced by fragmentation in high-energy heavv- 
ion reactions lead to the discovery of halo nuclei, such as 11 Li, about 20 years ago 
Nowadays a huge number of /3-unstable nuclei far from stability are being studied thanks 
to further technical improvements. Unstable nuclei far from stability are known to play an 
important role in nucleosynthesis. Detailed studies of the structure and their reactions will 
have unprecedented impact on astrophysics 0. 

The first experiments with unstable nuclear beams aimed at determining nuclear sizes 
by measuring the interaction cross section in high energy collisions Successive use of 
this technique has yielded nuclear size data over a wide range of isotopes. Other techniques, 
e.g. isotope-shift measurements, have allowed to extract the charge size. The growth of a 
neutron skin with the neutron number in several isotopes have been deduced from nuclear- 
and charge- size data 0]. 
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Undoubtedly, electron scattering off nuclei would provide the most direct information 
about charge distribution, which is closely related to the spatial distribution of protons 
0). A technical proposal for an electron-heavy-ion collider has been incorporated in the 



GSI/Germany physics program |5j]. A similar program exists for the RIKEN/ Japan facility 

jfj. In both cases the main purpose is to study the structure of nuclei far from the stability 
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line. The advantages of using electrons in the investigation of the nuclear structure are 
mainly related to the fact that the electron-nucleus interaction is relatively weak. For 
this reason multiple scattering effects are usually neglected and the scattering process is 
described in terms of perturbation theory. Since the reaction mechanism in perturbation 
theory is well under control the connection between the cross section and quantities such as 
charge distributions, transition densities, response functions etc., is well understood 0|. 

Under the impulse approximation, or plane wave Born approximation, the charge form 
factor can be determined from the differential cross section of elastic electron scattering. 
Since the charge distribution, p c h(r), is obtained from the charge form factor by a Fourier 
transformation, one can experimentally determine p c h{^) by differential cross-section mea- 
surements covering a wide range of momentum transfer q. One can obtain information on 
the size and diffuseness when the charge form factor is measured at least up to the first 
maximum. To do this within a reasonable measuring time of one week, a luminosity larger 
than 10 26 cm _2 s _1 is required, for example for the 132 Sn isotope 

On the theoretical side the difference between the proton and neutron distributions can 
be obtained in the framework of Hartree-Fock (HF)method (see for example [H) or Hartree- 
Fock-Bogoliubov (HFB) method (see for example jHEJ). As a rule of thumb, a theoretical 
calculation of the nuclear density is considered good when it reproduces the data on elastic 
electron scattering. But some details of the theoretical densities might not be accessible 
in the experiments, due to poor resolution or limited experimental reach of the momentum 
transfer q. 
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II. ELASTIC ELECTRON SCATTERING 



In the plane wave Born approximation (PWBA) the relation between the charge density 
and the cross section is given by 



ffi/pwBA ~ 1 + (2E/M A ) sin 2 (9/2) 1 ch W 1 ' U 

where a a/ = (Z 2 a 2 /4E 2 ) cos 2 (0/2) sin" 4 (0/2) is the Mott cross section, the term in the 
denominator is a recoil correction, E is the electron total energy, Ma is the mass of the 
nucleus and 9 is the scattering angle. 

The charge form factor F c h (q) for a spherical mass distribution is given by 

POO 

F c h (q) = 47r / dr r 2 j (qr) p ch (r) , (2) 
Jo 

where q = 2k sin (9/2) is the momentum transfer, hk is the electron momentum, and 
E = ^Jh 2 k 2 c 2 + ffigC 4 . The low momentum expansion of eq. El yields the leading terms 
F c h (q) /Z = 1 — q 2 (r 2 h ) /6 + • • • . Thus, a measurement at low momentum transfer yields a 
direct assessment of the mean square radius of the charge distribution, (r 2 h ) 1 ^ 2 . However, as 
more details of the charge distribution is probed more terms of this series are needed and, 
for a precise description of it, the form factor dependence for large momenta q is needed. 

A theoretical calculation of the charge density entering eq. 121 can be obtained in many 
ways. Let p p (r) and p n (r) denote the point distributions of the protons and the neutrons, 
respectively, as calculated, e.g. from single-particle wavefunctions obtained from an average 
one-body potential well, the latter in general being different for protons and neutrons. If 
Je p (r) and /# n (r) are the spatial charge distributions of the proton and the neutron in the 
non-relativistic limit, then the charge distribution of the nucleus is given by 

Pch (r) = J p p (r') f Ep (r - r') dV + J p n (r') f En (r - r') dV. (3) 

The second term on the right-hand side of eq. H3 plays an important role in the inter- 
pretation of the charge distribution of some nuclear isotopes. For example, the half-density 
charge radius increases 2% from 40 Ca to 48 Ca, whereas the surface thickness decreases by 
10% with the result that there is more charge in the surface region of 40 Ca than of 48 Ca 

EH- 

This also results that the rms charge radius of 48 Ca is slightly smaller than that of 
40 Ca. The reason for this anomaly is that the added fa/2 neutrons contribute negatively to 
the charge distribution in the surface and more than compensate for the increase in the rms 
radius of the proton distribution. 

For the proton the charge density Je p (r) in eq. El is taken as an exponential function, 
what corresponds to a form factor GE P (q) = fEp(q) — (1 + q 2 /A 2 )~ 1 . For the neutron a 
good parametrization is GE n (q) = - V^njG E P (q) / '(1 +pr), where p n is the neutron magnetic 
dipole moment and r = q 2 /4m jy. One can use A 2 = 0.71 fm -2 (corresponding to a proton 
rms radius of 0.87 fm) and p = 5.6, which Galster et al. |12] have shown to reproduce 
electron-nucleon scattering data. 

Eqs. ^Olare based on the first Born approximation. They give good results for light nuclei 
(e.g. 12 C) and high-energy electrons. For large-Z nuclei the agreement with experiments 
is only of a qualitative nature. The effects of distortion of the electron waves have been 
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studied by many authors (see, e.g. ref. 0,0, 15l|). More important than the change in the 
normalization of the cross section is the displacement of the minima. It is well known that a 
very simple modification can be done in the PWBA equation that reproduces the shift of the 
minima to lower g's. One replaces the momentum transfer q in the form factor of eq. ^with 
the effective momentum transfer q e ff = q (1 + 3Ze 2 /2R c hE), where E is the electron energy 
and R ~ 1.2 A -1 / 3 fm. This is because a measurement at momentum transfer q probes in 
fact |F(g)| 2 at q = q e ff due to the attraction the electrons feel by the positive charge of 
the nucleus. This expression for q e ff assumes a homogeneous distribution of charge within 
a sphere of radius R. 

A realistic description of the elastic electron scattering cross section requires the full 
solution of the Dirac equation. The Dirac equation for the elastic scattering from a charge 
distribution can be found in standard textbooks, e.g. [3]. Numerous DWBA codes based 
on Dirac distorted waves have been developed and are public. 



III. SKINS AND HALOS 
A. Neutron Skins 

Appreciable differences between neutron and proton radii are expected |l7[ to characterize 
the nuclei at the border of the stability line. The liquid drop formula expresses the binding 
energy of a nucleus with N neutrons and Z protons as a sum of bulk, surface, symmetry and 
Coulomb energies E/A = -a v A + a s A 2 /' 3 + S(N - Zf/A + a c Z 2 /A 1 /' 3 ± a p A _1 / 2 ,where a v , 
as, a p , S and ac are parameters fitted to the experimental data of binding energy of nuclei. 
This equation does not distinguish between surface (S) and volume (V) symmetry energies. 
As shown in ref. , this can be achieved by partitioning the particle asymmetry as N — Z = 
N s — Z s + N v — Zy. The total symmetry energy S then takes on the form 5* = Sy(N v — 
Zy ) 2 1 A + Ss ( Ns — Z.q) 2 1 A 2 ! 3 . Minimizing under fixed N — Z leads to an improved liquid drop 
formula 18] with the term S V (N - Z) 2 /A replaced by S V {N - Z) 2 /A[l + (S v /S s ) A' 1 / 3 ). 
The same approach also yields a relation between the neutron skin R np = R n — R p , and 
S s , S v , namely 18\R n - R p /R = (A/6NZ)[N -Z- (a c /12S v ) ZA 2 / 3 ]/[l + (S s /S v ) A 1 / 3 ], 
where R = (R n + R p )/2. 

Here the Coulomb contribution is essential; e.g. for N = Z the neutron skin R np is 
negative due to the Coulomb repulsion of the protons. A wide variation of values of Sy and 
Ss can be found in the literature. These values have been obtained by comparing the above 
predictions for energy and neutron skin to theoretical calculations of nuclear densities and 
experimental data on other observables [Til . 19l |. For heavy nuclei, with A>1, NZ ~ A 2 /2, 



and using a c = 0.69 MeV, S v = 28 MeV and R = 1.2A 1 ' 3 fm, one gets 

R np = R n -R p ~ 0.4 f^j (5 - 2.05 x Kr 3 ZA~ 1/3 ) fm, (4) 

where 5 = (N — Z) j A is the asymmetry parameter. If one assumes that the central densities 
for neutrons and protons are roughly the same and that they are both described by a uniform 
distribution with sharp-cutoff radii, R n and R P) one finds R np ~ 0.8A 1,/3 5 fm. Of course, the 
sharp sphere model is too simple. 

On the experimental front, a study of antiprotonic atoms published in reference |2o[ 
obtained the following fitted formula for the neutron skin of stable nuclei in terms of the 
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root mean square (rms) radii of protons and neutrons 

Ar np = (r 2 n ) 1/2 - (r 2 p ) 1/2 = (-0.04 ± 0.03) + (1.01 ± 0.15) 5 fm. (5) 

The relation of the mean square radii with the half-density radii is given by (r 2 ) = 3R 2 /5 + 
77r 2 a^/5, where a n is the diffuseness parameter. For heavy nuclei, assuming a n = a p ^ 
R n , R p , one gets the same linear dependence on the asymmetry parameter as in eq. |H 

Eq. |S]can be used as the starting point for accessing the dependence of electron scattering 
on the neutron skin of heavy nuclei. Applying it to calcium isotopes as an example, one 
obtains that the neutron skin varies from —0.15 fm for 35 Ca (proton-rich with negative 
neutron skin) to 0.25 fm for 53 Ca. A negative neutron skin obviously means an excess of 
protons at the surface. 

For heavy nuclei the charge and neutron distributions can be described by a Fermi dis- 
tribution. The diffuseness is usually much smaller than the half-density radius, a Pjn -C R p , n - 
The neutron skin is then given by R np = R n — R p ~ a/5/3 Ar np . One can assume further 
that the nuclear charge radius is given by R p = 1.2 A 1 ' 3 fm - R np /2 with Ar np given by eq. 

The first and second minima of the form factors occur at q\ = 4.49/ R p and q 2 = 7.73/ R p , 
respectively, corresponding to the zeroes of the transcendental equation tan (qR p ) = qR p . 

The linear dependence of R p with the neutron skin (and with the asymmetry parameter 
5), also imply a linear dependence of the position of the minima, 

-l 

fm -1 , and q 2 = 1.72 q 1 . (6) 

For 100 Sn the first minimum is expected to occur at q 1 = 0.806 fm" 1 = 159 MeV/c, while 
for 132 Sn it occurs at q x = 0.754 fm -1 = 149 MeV/c. 

The variation of qi with the neutron skin of neighboring isotopes, Ag x ~ 2/A 8 ^ 3 fm -1 , is 
too small to be measured accurately. The first minimum, q±, changes from 220 MeV/c for 
35 Ca to 204 MeV/c for 53 Ca, an approximate 7%, which is certainly within the experimental 
resolution. Of course, sudden changes of the neutron skin with 5 might happen due to shell 
closures, pairing, and other microscopic effects. 

To be more specific, let us assume that a reasonable goal is to obtain accurate results for 
the charge radius (r 2 ) 1 ^ 2 , so that 5 (r 2 ) 1 ^ 2 < 0.05 fm. This implies that the measurement of 
qi has to be such that (Aq 1 /q 1 ) < qi [fm -1 ] %, with qi in units of fm -1 and the right-hand 
side of the inequality yielding the percent value. For 53 Ca, qi = 1.11 fm -1 meaning that the 

1 /2 

experimental resolution on the value of qi has to be within 1% if 5 (r 2 ) < 0.05 fm is a 
required precision. Of course, the ultimate test of a given theoretical model will be a good 
reproduction of the measured data, below and beyond the first minimum. 

B. Neutron halos 

Elastic electron scattering will be very important to determine charge distributions in 
proton-rich nuclei. This will complement the basic information on the charge distribution 
in, e.g. 8 B, obtained in nucleon knockout reactions 2l]]. For light nuclei composed by a 
core nucleus and an extended distribution of halo nucleons, the nuclear matter form factor 
can be fitted with the simple expression F(q) = (1 — g) exp (— q 2 a\/4) + g/(l + a^q 2 ), with 
the density normalized to one. g is the fraction of nucleons in the halo. In this expression 
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the first term follows from the assumption that the core is described by a Gaussian and 
the halo nucleons by an Yukawa distribution. Taking n Li as an example, the following 
set of parameters can be used g = 0.18, a± = 2.0 fm and ci2 = 6.5 fm. Although only few 
nucleons are in the halo they change dramatically the appearance of the squared form factor. 
Even when the individual contribution of the halo nucleons is small and barely visible in a 
linear plot of the matter distribution, it is very important for the form factor of the total 
matter distribution. It is responsible for the narrow peak which develops at low momentum 
transfers. This signature of the halo was indeed the motivation for the early experiments 
with radioactive beams. The narrow peak was observed in momentum distributions following 
knockout reactions 

Elastic electron scattering will not be sensitive to the narrow peak of |F(g)| 2 at small 
momentum q, but to the form factor of the charge distribution, |F c ^(g)| 2 . The determination 
of this form factor will tell us if the core has been appreciably modified due to the presence 
of the halo nucleons. 

In order to explain the spin, parities, separation energies and size of exotic nuclei con- 
sistently a microscopic calculation is needed. One possibility is to resort to a Hartree-Fock 
(HF) calculation. Unfortunately, the HF theory cannot provide the predictions for the sep- 
aration energies within the required accuracy of hundred keV. I have used a simple and 
tractable HF method to generate synthetic data for the charge-distribution of n Li. 
Details of this method is described in ref. j2^|. Assuming spherical symmetry, the equation 
for the Skyrme interaction can be written as 

h 2 

-V + V(r) 



2m* (r) 



^(r)=e^(r) (7) 



where m*{r) is the effective mass. The potential V(r) has a central, a spin-orbit and a 
Coulomb term. The central potential is multiplied by a constant factor / only for the last 
neutron configuration: 

t/ t \ f\r ( \ / / 7^ 1 f° r l as t neutron configuration , . 

Vcentrai(r) _ }V HF {r), \ f = 1 otherwise. (8) 

As the effective interaction, a parameter set of the density dependent Skyrme force, so 
called BKN interaction 24]|, is adopted. The parameter set of BKN interaction has the 
effective mass m*/m =1 and gives realistic single particle energies near the Fermi surface 
in light nuclei. The original BKN force has no spin-orbit interaction. In the calculations, 
the spin-orbit term was introduced in the interaction so that the single-particle energy of 
the last neutron orbit becomes close to the experimental separation energy. In this way, the 
asymptotic form of the loosely-bound wave function becomes realistic in the neutron-rich 
nucleus. In order to obtain the nuclear sizes, the rms radii of the occupied nucleon orbits 
are multiplied by the shell model occupation probabilities, which are also obtained in the 
calculations. 

The elastic form factor for the matter distribution obtained in the HF calculations are 
very close to the ones calculated by the empirical formula given at the beginning of this 
section. The lack of minima, and of secondary peaks (as in the empirical formula), make 
it difficult to extract from \F c h(q)\ 2 more detailed information on the charge- density profile. 
For example, in the case of 6 Li a good fit to experimental data was obtained with [25| 
\F ch (q)\ 2 oc exp(-aV) - Cq 2 exp(-6 2 g 2 ), with a = 0.933 fm, b = 1.3 fm, and C = 0.205. 
However, the data cannot be fitted by using a model in which the nucleons move in a 
single-particle potential. 
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C. Proton Halos 



I will consider 8 B as a prototype of proton halo nucleus. This nucleus is perhaps the most 
likely candidate for having a proton halo structure, as its last proton has a binding energy of 
only 137 keV. The charge density for this nucleus can be calculated in the framework of the 
Skyrme HF model. I use here the results obtained in ref. [26j], where axially symmetric HF 
equations were used with SLy4 (22J Skyrme interaction which has been constructed by fitting 
the experimental data on radii and binding energies of symmetric and neutron-rich nuclei. 
Pairing correlations among nucleons have been treated within the BCS pairing method. The 
form factor squared for the charge density in 8 B is calculated. 

The width of the charge form factor squared is A c h = 0.505 fm _1 = 99.6 MeV/c. The 
corresponding values for the neutron and the total matter distributions are, respectively, 
A n = 0.512 fm" 1 = 101 MeV/c and A tot = 0.545 fin" 1 = 108 MeV/c. This amounts to 
approximately 10% differences of matter and charge form factors in 8 B. 

The proton halo in 8 B is mainly due to the unpaired proton in the P3/2 orbit. It is clear 
that for such a narrow halo the size of the nucleon also matters. A slice of the nucleon is 
included in a thin spherical shell of radius r and thickness dr from the center of the nucleus. 
If the position of the nucleon is given by R, the part of the proton charge included in the 
spherical shell is given by 

dp c h = 2nr 2 dr / d6 p p (x) sin 9, (9) 
Jo 

where p p (x) is the charge distribution inside a proton at a distance x from its center. 
The coordinates are related by x 2 = r 2 + R 2 — 2rRcos8. The contribution to the nuclear 
charge distribution from a single-proton in this spherical shell is thus given by J\f p (R, r) = 
dp ch /A-nr 2 dr. 

Assuming that the charge distribution of the proton is described either by a Gaussian or 
a Yukawa form, the integral in eq. |H1 can be performed analytically, yielding 

for a Gaussian dist. 
for a Yukawa dist., 

(10) 

where a is the proton radius parameter. 

The charge distribution at the surface of a heavy proton- rich nucleus, 5p c h (r), may be 
described as a pile-up of protons forming a skin. Let Hi be the number of protons in the 
skin and Ri their distance to the center of the nucleus. One gets 

Sp ch {r) = J2 n iK(^ R i)- (11) 

i 

Assuming Ri to be constant, equal to the nuclear charge radius R, and using eq. [TQ] it 
is evident that while the density at the surface increases, its size R and width a, remain 
unaltered. The form factor associated with this charge distribution is given by 

._. t N 47r ^-^ [°° . . , 2N ^-^ sin (qRi) _ x 
0F ch (q) = — y rti / dr r N p (r, Ri) sin (qr) = exp l-qa 1 > m , (12) 

Q i Jo i Q R i 
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where the last result is for the Gaussian distribution. An analytical expression can also be 
obtained for the Yukawa distribution. For Ri = R, expression ^] shows that the increase of 
density in the skin does not change the shape of the form factor, or of the cross section, but 
just its normalization. The fall down of the cross section is determined by a alone, and not 
by rii. If the charge of additional protons is distributed homogeneously across the nucleus 
including the skin, the normalized form factor will not change, except for a small change in 
R. 

For a proton halo nucleus it is more appropriate to replace Yli n % ~~ * 
Ati J dR R 2 Af p (r, R) 5p c h (R), where 5p c h (R) is the density change created by the extended 
wavefunction of the halo protons. One then recast eq. ED in the form 

Ait r°° 

6F ch (q) = — exp (-qa 2 ) / dR R Sp ch (R) sin (qR) . (13) 
Q Jo 

The shape of the form factor is here dependent not only on the proton size but also on the 
details of the halo density distribution. For 8 B, the halo size is determined by the valence 
proton in a P3/2 orbit. The density 5p c h (R) due this proton can be calculated with a Woods- 
Saxon model. Using the same potential parameters as in ref. [2^ I compare the form factor 
\5F ch (q)\ 2 to the charge form factor \F ch (q)\ 2 . The halo is found to contribute to a narrow 
form factor. However, in contrast to the neutron halo case, the charge form factor of 8 B 
does not show a pronounced influence of the halo charge distribution. 

The rms radius of the charge distribution can be calculated from (r 2 h ) = 

—6 [dF c h/d (q 2 )} 2=0 . Applying this relationship to the 8 B charge form factor we get (r 2 ^ 1 ^ 2 = 

2.82 fm which is close to the experimental value (r 2 ^]/ 2 = 2.76 ± 0.08 fm. The shape of the 
charge form factor can also be described by a Gaussian distribution with radius parameter 
a = 2.30 fm. In contrast to the case of n Li, the proton halo in 8 B does not seem to build 
up a two-Gaussian shaped form factor. This observation also seems to be compatible with 
the momentum distributions of 7 Be fragments in knockout reactions using 8 B projectiles in 
high energy collisions j2l|. Electron scattering experiments will help to further elucidate 
this property of proton halos. 



IV. INVERSE SCATTERING PROBLEM 

In PWBA the inverse scattering problem can be easily solved. It is possible to extract 
the form factor from the cross section and then, with an inversion of the Fourier transform, 
to get the charge density distribution 

1 f°° 

Pch (r) = — — / dq q 2 j (qr) F ch (q) . (14) 
^ Jo 

The PWBA approximation can be justified only for light nuclei (e.g. 12 C) in the region far 
from the diffraction zeros. For higher Z values the agreement with experiment is only of a 
qualitative nature. 

It is very common in the literature to use a theoretical model for p c h (r), e.g. the HF 
calculations discussed in the previous sections and compare the calculated F(q) with ex- 
perimental data. When the fit is "reasonable" (usually guided by the eye) the model is 
considered a good one. However, whereas the theoretical p c h (r) can contain useful informa- 
tion about the central part of the density (e.g. bubble-like nuclei, with a depressed central 



7 



density), an excellent fit to the available experimental data does not necessarily mean that 
the data is sensitive to those details. The obvious reason is that short distances are probed 
by larger values of q. Experimental data from electron-ion colliders will suffer from limited 
accuracy at large values of q, possibly beyond q = 1 fm _1 . 

In order to obtain an unbiased "experimental" p c h (r) one usually assumes that the density 
is expanded as p c h (r) = J2^Li a nfn ( r ) ? where the basis functions f n ( r ) ar e drawn from any 
convenient complete set and the expansion coefficients a n are adjusted to reproduce the 
differential elastic cross section. The corresponding Fourier transform then takes the form 

00 /»oo 

P{q) = F ch (q) = a ^ n (?) ' f n (?) = 4vr / dr r 2 j (qr) f n (r) . (15) 

n=l J 

Evidently the sum in n has to be truncated and this produces an error in the determination 
of the charge density distribution. Another problem is that, as shown by eq. the solution 
of the inverse scattering problem requires an accurate determination of the cross section 
up to large momentum transfers. Electron scattering experiments in electron-ion colliders 
will be performed within a limited range of q and this will produce an uncertainty in the 
determination of the charge density distribution. 

Two bases have been found useful j2^| in the analysis of electron or proton scattering 
data. The present discussion is limited to spherical nuclei, but generalizations to deformed 
nuclei can be done. The Fourier-Bessel (FB) expansion (i.e. with /„ taken as spherical 
Bessel functions) is useful because of the orthogonality relation between spherical Bessel 
functions 

rRmax 1 

/ dr r 2 ji (q n r) j, (q m r) = -R 3 C \j t+1 (q n R mSLX )} 2 8 nm , (16) 
Jo * 

where the q n are defined such as ji (g„-R max ) = 0. The FB basis implies that the charge 
density p c h(r) should be zero for values of r larger than -R max . For example, the basis can 
be defined as follows 

fn (r) = jo (q n r) 6 (i? max - r) , f n (q) = 4?r ( ^ ^ max j Q ( q R max ) , (17) 

q -q n 

where 9 is the step function, -R max is the expansion radius and q n = nii/R max . 

In principle it is possible to obtain the a n coefficients measuring directly the cross section 
at the q n momentum transfer. If the form factor (J2J) is known at q n , the coefficients a n can 
be obtained inserting (|16|) and (|17j) in the definition (J2J) of the form factor, leading to 

n _ RciMn) _ (lg) 



[Ji (^n^max)]' 

In general the cross sections are measured at q values different from q n . Using the 
expansion (j!7)l of the charge density one finds for the form factor the relation 

Att ( l) n 

Fch (q) = — V" a n ^ 5- sin (qR max ) . (19) 

q V q q n 



By fitting the experimental F c h{q) one obtains the a n parameters and reconstruct the nuclear 
charge density. Not all a n 's are needed. Since the integral of the density, or F (q = 0), is 
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fixed to the charge number there is one less degree of freedom. Also, densities tend to zero 
at large r. Thus another condition can be used, e.g. that the derivative of the density is 
zero at R maiX . Thus, when one talks about n expansion coefficients one means in fact that 
only n — 2 coefficients need to be used in eq. For experiments performed up to g max the 
number of expansion coefficients needed for the fit is determined by n max ~ ^max-Rmax/ 71 "- 

A disadvantage of the FB expansion is that a relatively large number of terms is often 
needed to accurately represent a typical confined density, e.g. for light nuclei. One can use 
other expansion functions which are invoke less number of expansion parameters, e.g. the 
Laguerre- Gauss (LG) expansion, 



where a = r//3, 7 = q(3/2, and L v n is the generalized Laguerre polynomial. Another pos- 
sibility is to use an expansion on Hermite (H) polynomials. In both cases, the number of 
terms needed to provide a reasonable approximation to the density can be minimized by 
choosing (3 in accordance with the natural radial scale. For light nuclei /3 = 1 — 2 fm can 
be chosen, consistent with the parametrization of their densities. Then the magnitude of a n 
decreases rapidly with n, but the quality of the fit and the shape of the density are actually 
independent of (3 over a wide range. 

For real data, the expansion coefficients obtained by minimizing 



where y(qi, a n ) is the fitted value of the cross section (form factor) with a set of coefficients 
a n and y~i are the experimental data at momentum qi with uncertainty o",. 

Increasing the number of coefficients does not improve the quality of the fit. It only 
produces more oscillations of the density. The reason is that terms with larger n's are only 
needed to reproduce the data at larger values of momentum transfer. The fit to the data 
for q < g max is not affected but the presence of these new terms introduces oscillations in 
the charge distribution. A possible fix to this problem is to include pseudodata in addition 
to experimental data. This method is well known in the literature j2^|. The pseudodata 
are used to enforce constraints and to estimate the incompleteness error associated with the 
limitation of experimental data to a finite range of momentum transfer. 

The author is grateful to Haik Simon for beneficial discussions. This work was supported 
by the U. S. Department of Energy under grant No. DE-FG02-04ER41338. 
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